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Abstract 

^ , The annual occurrence of many infectious diseases remains a constant burden to 

a^ _ 

' public health systems. The seasonal patterns in respiratory disease incidence observed 

' in temperate regions have been attributed to the impact of environmental conditions 

■ on pathogen survival. A model describing the transmission of an infectious disease by 

means of a pathogenic state capable of surviving in an environmental reservoir out- 
. side of its host organism is presented in this paper. The ratio of pathogen lifespan 

I to the duration of the infectious disease state is found to be a critical parameter in 

determining disease dynamics. The introduction of a seasonally forced pathogen inac- 
tivation rate identifies a time delay between peak pathogen survival and peak disease 
incidence. The delay is dependent on specific disease parameters and, for influenza, 
decreases with increasing reproduction number. The observed seasonal oscillations are 
found to have a period identical to that of the seasonally forced inactivation rate and 
which is independent of the duration of infection acquired immunity. 
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1 Introduction 



Many diseases exhibit seasonal cycles in incidence data, most notably childhood diseases 
such as measles and rubella (Anderson and May, 1991, Dowcll, 2001). For such diseases, 
it is widely accepted that the observed cycles are linked to the timing of school terms and 
the subsequent increase in contact rates among the immunologically naive child population 
(Keeling and Rohani, 2001, Stone et al. 2007). The mechanism driving the seasonal occur- 
rence of influenza and other respiratory infections is less well understood. While increased 
contact rates due to indoor crowding during winter undoubtedly facilitate greater disease 
transmission, it is more likely a contributing factor and not the driving mechanism (Lofgren 
et al. 2007). Seasonal influenza places a considerable burden on public health systems with 
annual global incidence in the range 5 - 15% of the population, resulting in up to 500,000 
deaths (Stohr, 2002). Targeted interventions to reduce this burden could be significantly 
improved if the seasonal stimulus was better understood. 

The cause of seasonality must be attributed to annual changes associated with either 
the host or the infectious disease pathogen (Crassly and Fraser, 2006). In addition to in- 
creased contact rates during winter, higher transmission rates have been linked to changes 
in the human immune system. Possible variations in human immune function have been 
attributed to fluctuations in melatonin secretion regulated by the annual light-dark cycles 
and to deficiencies in vitamin D during the winter months (Cannell et al. 2008, Lofgren 
et al. 2007, Nelson and Demas, 1996). Experiments on mice in an environmentally con- 
trolled environment found that mice were more susceptible to infection during the winter 
months (Schulman and Kilbourne, 1963). For diseases transmitted through indirect path- 
ways (e.g. respiratory droplets, fomites, fungal spores, waterborne pathogens) the ability 
of the pathogen to survive outside of its host must play a vital role in the transmission 
process (Crassly and Fraser, 2006). This survival potential would be greatly infiuenced by 
environmental factors. Early studies on the survival of the infiuenza virus in air indicated 
that pathogen survival peaks at low relative humidity and low temperatures resulting in 
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increased viral transmission (Harper, 1961, Hemmes et al. 1960). Experiments to study the 
airborne transmission of influenza among Guinea pigs further supported this theory indicat- 
ing that transmission among subjects peaked in cold dry air (Lowen et al. 2007). However, 
while indoor relative humidity is minimized in winter, outdoor humidity peaks. More re- 
cent studies identify absolute humidity as a more likely seasonal driver, as both indoor and 
outdoor absolute humidity display cyclic behavior in temperate regions that minimizes in 
winter (Shaman and Kohn, 2009). Influenza virus survival increases with low absolute hu- 
midity leading to increased transmission in winter (Shaman et al. 2010). Other theories for 
increased viral inactivation in summer consider ultraviolet radiation and air toxicity, how- 
ever, little work has been undertaken to support such theories (Weber and Stilianakis, 2008). 

A crucial factor to consider when assessing the impact of environmental factors on sea- 
sonal transmission is the mode of disease transmission. For respiratory diseases three modes 
of pathogen transmission have been identified: droplet, contact and airborne transmission 
(Weber and Stilianakis, 2008). Droplet transmission occurs when, following an expiratory 
event by an infected individual, large pathogen-carrying droplets (diameter > 10 /im) are 
deposited directly onto the mucous membranes of a susceptible person. The direct nature 
of such a transmission pathway precludes a significant impact by environmental factors on 
either the droplet itself or its pathogen load. Conversely, the contact (through fomites or di- 
rect human-to-human contact) and airborne (small aerosol droplets with diameter < 10 /im) 
modes of transmission render the pathogen vulnerable to environmental conditions during 
prolonged periods spent external to its host organism. For the specific case of the infiuenza 
virus, the efficiency of contact transmission is determined by the survival rate of the pathogen 
on solid surfaces and human skin. The high inactivation rates observed on hands can limit 
the occurrence of transmission via the contact route, however, a continuous supply of fomites 
from infected individuals could possibly counteract this (Weber and Stilianakis, 2008). This 
mode of transmission is greatly infiuenced by human behavior (e.g. cough etiquette, hand 
washing) and thus a seasonal variation in its impact would be difficult to quantify. Airborne 
droplets remain exposed to the ambient environmental conditions for prolonged periods of 

3 



time and could act as an important agent driving seasonal disease incidence. For respiratory 
diseases mediated by airborne droplets and for a spatially homogeneous distribution of sus- 
ceptibles, pathogen removal is the result of three distinct processes: gravitational settling, 
pathogen inactivation and inhalation (Stilianakis and Drossinos, 2010). Environmental fac- 
tors are capable of influencing the first two of these processes (Shaman and Kohn, 2009). 
The gravitational settling rate of a droplet is determined from its diameter, which is influ- 
enced by evaporation effects determined by the ambient air properties. However, a review of 
experimental studies by Shaman and Kohn (2009) found that there is insufficient evidence 
to support such a hypothesis. They concluded that seasonal variations in environmental 
conditions are most likely incorporated through the pathogen inactivation term. 

The seasonality of infectious diseases is usually incorporated in deterministic models 
through a time dependent transmission rate. Outbreaks are typically simulated using a si- 
nusoidal function or, in the case of childhood diseases such as measles, a step function is 
employed to represent school terms (Crassly and Fraser, 2006, Keeling and Rohani, 2001). 
Models incorporating natural birth and death processes display damped oscillations towards 
an endemic disease state (Fisman, 2007). However, if the period of the applied seasonal 
forcing is close to the intrinsic period of the damped oscillations then the two effects can res- 
onate to produce large amplitude seasonal oscillations (Dushoff et al. 2004). An important 
consideration in any analysis of seasonal influenza is the process of antigenic drift, whereby 
continuous small changes in the virus requires the production of a new vaccine each year. 
This phenomenon can be integrated into the standard deterministic models by allowing re- 
covered individuals to lose their acquired immunity after a specifled period of time (Dushoff 
et al. 2004). An important factor neglected in the standard deterministic models is the 
mechanism driving the seasonal variation. If the theory that the impact of annual variations 
in environmental conditions on pathogen survival/inactivation is indeed valid then models 
should incorporate this phenomenon. In this paper, we present such a model. 

In Section [2] we analyze a model for the spread of an infectious disease by means of an 
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intermediate free-living pathogenic state, which is exposed to the ambient environmental 
conditions. This general model, applicable to a variety of free-living organisms (e.g. viruses, 
bacteria, fungi, protozoa), was introduced by Anderson and May (1981) to describe indirect 
disease transmission between free-living microparasites and their invertebrate hosts. Varia- 
tions of the basic model have been employed to describe the spread of a pathogen through a 
generalized environmental state (Li et al. 2009), the spread of a waterborne bacteria (Tien 
and Earn, 2010), bacterial and prion disease in livestock and wildlife (Nieuwhof 2009, Miller 
2006), and the transmission of respiratory diseases by airborne droplets (Stilianakis and 
Drossinos, 2010). The model has even been adapted to describe the spread of fungal spores 
in a vineyard (Burie, 2006) and the point release of an infectious agent (Reluga, 2004). With 
such a wide range of applications and descriptive abilities the intrinsic characteristics of the 
infectious agent are a vital model component and can radically impact the dynamics. A 
pathogen characteristic of primary importance is its ability to survive for prolonged periods 
outside of its host, which is directly influenced by environmental conditions and, thus, can- 
not be ignored in the quest to identify seasonal disease drivers. Therefore, we first analyse 
the general model and consider its behavior for both short-lived and long-lived pathogenic 
states. An important outcome of this analysis is that, in the case of short-lived pathogens, 
a quasi-steady state exists whereby the pathogen dynamics can be described in terms of the 
infected population alone. This quasi-steady approximation can also be applied to an au- 
tonomous model describing seasonal outbreaks. This model is analysed in Section [3l where 
we consider the consequences of a seasonal variation in the pathogen inactivation rate and 
how it impacts disease incidence. 

2 A model for indirect transmission of an infectious 
disease 

Many infectious diseases are primarily transmitted by means of an intermediate environ- 
mental reservoir, in which free-living pathogens are capable of surviving outside of the host 
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organism. For example, waterborne outbreaks can persist through the shedding of pathogens 
by an infected individual into a water source which is then ingested by a susceptible individ- 
ual. Similarly, respiratory diseases can be transmitted from person to person via pathogen 
loaded airborne droplets or fomites expelled into the environment by infected individuals 
during expiratory events (e.g. coughing or sneezing). The efficiency of such transmission 
pathways will depend on the ability of the pathogen to survive in the intermediate reservoir 
(e.g. air, water). In this section we present a generalized model for the transmission of an 
infection by means of such a pathogen reservoir. 

Consider a closed population of individuals, of which S{t) are susceptible, I{t) are 
infected and R(t) are recovered. Infected individuals shed pathogens into an intermediate 
reservoir, which are free-living (outside of the host organism) and can transmit infection 
through contact with the susceptible population. Figure [H The total number of pathogens 
in the reservoir is P{t). A deterministic model for the disease dynamics is 



where the total human population A^ = S' + / + -Ris constant. The duration of time spent in 
the infected state is T/ = l/(/i + rf), where is the average infectious disease period and 
l/c? is a typical human lifespan. Infected individuals generate pathogens, with a lifespan of 
Tp = l/a, at a rate k. Infection acquired immunity is lost after a time period of l/a, with 
cr = if immunity is permanent. The transmission rate per pathogen is /3, determined from 
the contact rate of a susceptible with a pathogen and the probability of such a contact caus- 
ing infection. The equation describing pathogen dynamics assumes that pathogen numbers 
decrease exponentially for positive a. This could be modified to describe pathogens with 
more complicated life cycles. The pathogen removal rate a can be a complicated function of 
environmental factors, depending on the specific pathogen type and reservoir characteristics. 
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Figure 1: Infection pathway for the transmission of a disease by a free-hving pathogen P in 
a population of susceptible 5", infected / and recovered R individuals. 



but for now we assume it is constant. 



The general model presented here is apphcable to a wide variety of diseases such as 
cholera and influenza in humans (Tien and Earn 2010, Stilianakis and Drossinos 2010) and 
footrot and chronic wasting disease in animals (Nieuwhof 2009, Miller 2006). In the case 
of waterborne pathogens, this model has been expanded, specifically in the case of cholera, 
to include dose dependent infection rates (Codego 2001, Bertuzzo 2010) and a minimum 
infectious dose (Joh, 2009). To analyse the general framework of the basic model we omit 
such disease dependent complications, but will consider this issue in relation to influenza in 
Section 3.3. 

The model can be written in dimensionless form by scaling 

S,I^N, P^-N, t —„ (1) 

a IJ. + a 
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and the dimensionless system is 



^ = ^{l-S) + iil -S-I)- RiPS, (2) 
'^-RoPS-I, (3) 

where dimensionless numbers are defined as 

/3k d a n + d Tp 

Ro = —, — —17, <P = — — 1, W = — — 1, P= = 7^- 

a{iJ, + d) jji + d jjL + d a Ij 

The dimensionless number p represents the ratio of the pathogen lifespan to that of an 
infected individual. The dynamics of the pathogen, satisfying (jl]), will be determined by the 
magnitude of p. The magnitude of i?o represents the average number of secondary infections 
from a primary infection and is typically 0(1). Finally, and ip determine how quickly the 
pool of susceptible individuals is repopulated from natural birth or death processes and loss 
of acquired immunity respectively. Suitable initial conditions are 

S{0) = So, I{0) = l~So = Io, P{0) = Po, (5) 
where it is assumed that no pre-existing immunity is present in the population, -R(O) = 0. 

The system ©-(jl]) has two equilibrium states. A disease-free state Eq = (1, 0, 0) and an 
endemic state E^, where 

p_.c^ , P^-M (0 + ^)(i^o-l) i<P + ^){Ro-l) \ 

which exists when at least one of or t/^ is nonzero and Rq > 1. It can be easily shown 
by a linear analysis that the disease-free state is locally stable when Rq < 1 and the initial 
infection ([5]) dies out. Conversely, the endemic state is locally stable when i?o > 1 and the 
initial infection will spread approaching the equilibrium state E^ as t ^ oo. The equilibrium 
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point is a stable focus resulting in damped oscillations towards Ee- As we are interested in 
the long-term seasonal variation in the endemic state we will assume that Rq > 1 for our 
purposes. 



2.1 Short-lived pathogenic state 

If the duration of the infectious disease state is significantly longer than the lifespan of 
the pathogen, Tp ^ Tj, then p 1. This implies that equation (jll) rapidly achieves an 
equilibrium compared with the rest of the system. We can assume that the left-hand side of 
equation is approximately zero, = 0, and employ the approximation 



(7) 



The system then reduces to the standard SIRS model 

dS 



H'^-S)+^{l-S-I)-RoIS, (8) 
^ = {RoS - 1)/, (9) 



with initial conditions 



3(0) = So, /(0) = Jo, P(0)=/o. (10) 

The quasi-steady state approximation employed above does not imply that ^ = 0, since 
P will change in time according to ([7]). A closed form solution is not possible for the SIRS 
system. However, the long term behavior is characterized by two possible equilibrium states. 
A disease-free state Eq = (1, 0) and the endemic state Ee where 

p -rc n - M (0 + ^)(i?o-i) \ 
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Figure 2: Numerical solution of the SIRS system OH])-® with parameter values (j) = ip = 0.1, 
i?o = 2 and p = 0.001. Initial conditions are ^(0) = 0.9, /(O) = 1 - S{Q) and P(0) = 10"^. 
(a) Long-time behavior of the infected population (b) Boundary layer region showing the 
small-time behavior of the infected population and pathogen number. 

and Pe = le from ([7]), in agreement with (|6]). It is trivial to show by a linear analysis that, 
in the limit t — )■ oo, the SIRS system approaches this stable endemic state when Rq > 1, 
which is shown by a numerical solution in Figure |2]^a). From f lTO|) . the quasi-steady state 
assumption requires -P(O) = Jq and the original initial conditions ([5]) are not satisfied, except 
in the unique case where Pq = Iq. This inconsistency arises because the approximation ([7]) is 
only valid when ^ ~ 0{1) and breaks down if P changes rapidly at any point in the domain 
such that is no longer small. There is an initial boundary layer for times of order p in 
which the approximation is invalid and the original scalings ([1]) are not appropriate. Figure 
m^b). In this region ^ 1 and we rescale the time in equations ©-(jll) using t = pr to 
obtain equations for this boundary layer region 

p[<Pil-S)+^{l-S-I)-RoPS], 
p [RoPS - I] , 
I -P. 
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At leading order in p we then have 




dl 
dr 







dP 
dr 



I -P. 



which has solution 



S = S, 



I 



P 



Jo + (Po - /o)e 



The original initial condition P(0) = Pq is now satisfied and as r — oo we exit the boundary 
layer and we have P — ?■ Jq in agreement with the quasi-steady state approximation flTU]) . 
In theory, unknown constants generated from the solution of the SIRS model in the outer 
region would be fixed by matching the two solutions in the limits as t — )■ and r — )■ oo. In 
terms of the long-term behaviour of an infectious disease, we can conclude that the dynamics 
of the pathogen closely follow those of the infected population and a standard SIRS model 
is sufficient for a mathematical analysis. However, because the pathogen undergoes rapid 
changes, the initial transient interval may be important for short-term dynamics. 

2.2 Long-lived pathogenic state 

If the pathogen survives significantly longer than the duration of the infectious disease state, 
Tp ^ Tf, then p ^ 1. For convenience, in this section we set ip = so that immunity is 
permanent. This has the effect of slowing the rate at which the susceptible compartment is 
repopulated. However, the behavior of the solution is qualitatively similar and the simplified 
system serves to highlight the primary features. The endemic equilibrium of the full system 
©-(HI), with ip = 0, is now given by 



Ee — (5'e, le, Pe) 




(11) 
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Figure 3: Numerical solution of the system ([2])-(j4]), obtained with initial conditions >S'(0) = 
0.9, /(O) = 0.1 and P(0) = 0.01. Parameter values are = 0.1, = 0, Rq = 2 and p = 1000. 
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and the solution approaches i?e as t — )■ oo when Rq > 1. When p ^ 1 the right-hand side of 
equation (jl]) is approximately zero, (/ — P)/p ~ 0, and we can assume that P is in a steady 
state such that ^ = and P = Pq. The system reduces to 

^ = 0(1 - 5) - RoPoS, 
dt 

— = RqPqS — /, 
dt 

with initial conditions 

^(0) = ^0, /(O) = Jo. 
The linear system is easily solved to obtain 

S{t) = S* + {So - S*)e-^\ (12) 
Ht) = !* + - e-*) + (/o - ne-\ (13) 



where 



+ RqPq 



The solution of the reduced system clearly satisfies the initial conditions of the full system 
^ as t — > 0. However, as t — )• oo the solution does not approach the endemic equilibrium 
ffTTj) . At relatively large times, the solution approaches an alternative endemic state E* = 
{S*, I*, Pq), Figure [3l This inconsistency arises because, for times t > 0(p), P is no longer at 
a steady state and begins to grow towards its Ee value. Figure HI Under such circumstances 
the steady state assumption breaks down. To describe this large-time behaviour, we rescale 
the time t = p6 and the system becomes 

= ^{1 - S) - RoPoS, 

p dO 
p dO 

'-^^I-P 

de 
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Figure 4: Numerical solution of the system (l2])-(jl]), obtained with initial conditions 5'(0) = 
0.9, J(0) = 0.1 and P(0) = 0.01. Parameter values are = 0.1, ^ = 0, i?o = 2 and p = 1000. 
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Now, for p ^ 1 this can be approximated by 



which has solution 



0(1 -S)- RqPS = 0, 
RoPS -1 = 0, 

de 



Q^m T(n^ <PRoPiO) .... 



where P{0) satisfies the first-order equation 

dP (j)RoP{0) 



p{e). (16) 



de <p + RoPie) 

For this large time solution to match the inner one, (fT2|) - (fT3|) . we require 

s{e ^ 0) = s*, i{e ^ 0) = r, p{e ^ o) = Pq- 

which yields 

s{9 = 0) = -4--, i{e = 0) 



+ RoP^ ' ' <P + RoPo 

in agreement with ( IT4|) . Equation ( !T6|l cannot be solved explicitly for P{0), however, an 
implicit solution can be obtained 

Po ( <P + RoP-<PRo Y" ^ eiRo-i) 
P \(l) + RoPo-(l)RoJ 

Now, as 6' — )• oo we find that the limiting value of P{9), which we denote as Poo, niust satisfy 



+ PoPoo-0Po "^° 

) + Po^o - (pRo 



0, 



which yields 



0(Po - 1) p 

Kq 



15 



It follows immediately from (fTS!) that S ^ Se and / — i- Je as — )■ oo and the large time 
solution approaches the endemic equilibrium ( ITTj) as required, Figure HI 

2.3 Applications of the model 

The model can be used to describe the airborne spread of a respiratory disease by aerosol 
droplets. The total number of airborne pathogens is P{t) = NpD(t), where D{t) is the total 
number of droplets with active pathogens and Np is the number of pathogens per droplet at 
the time of expulsion (Stilianakis and Drossinos, 2010). The lifespan of droplets is identical 
to that of pathogens (Tp = To = l/a) as droplets and pathogens are removed by identical 
processes since they are intrinsically linked. Gravitational settling removes droplets from 
the environment, which invariably results in the removal of the droplet's pathogen load. 
Similarly, the inhalation of a droplet by a population member removes pathogens, however, 
this process had little impact on droplet numbers and can be neglected (Robinson et al. 
2012). In addition, a droplet is removed from the infection pathway through pathogen 
inactivation (i.e. an airborne droplet carrying inactivated pathogens cannot cause infection). 
The dimensional system describing the infection dynamics is 



where /3d = f3Np and = ^ are the droplet transmission and generation rates respectively. 
Scaling D ~ — A^, the dimensionless system is 



dS_ 

m 

(H 

It 
dD 

~dt 



d{N - S) + a{N - S - I) - ^DS, 
^DS-{f, + d)I, 



Kdl - aD, 



dS 



4>{1 -S) + ij{l-S-I)- RoDS, 



(17) 



dt 
dl 



RqDS - /, 



(18) 



dt 
dD 



I-D, 



(19) 
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where all other scales and dimensionless numbers are defined as before. For the particular 
case of influenza the infectious period is approximately 5 days, fi = 0.2/day, and an average 
human lifespan is 70 years, d — 4 x 10~^/day. Airborne droplets (of diameter 4 /im) are 
removed, through the processes of gravitational settling and pathogen inactivation, at a rate 
of a = 37.44/day (Robinson et al. 2012). The droplet generation rate Kd is based on the 
number of pathogen loaded droplets emitted during a cough. Generation rates per cough 
are taken as 160/day (Nicas et al. 2005). The daily generation rate is then obtained by con- 
sidering a 200-fold increase for a sneeze (Nicas et al. 2005) and a total of 11 sneezes and 360 
coughs per day (Atkinson and Wein, 2008), which yields — 4.1 x 10^/day. Loss of immu- 
nity is typically associated with the emergence of new viral strains, with previously infected 
individuals reverting to a susceptible state after approximately 5 years, cr = 5.5 x 10~^/day 
(Truscott et al. 2012). 

For an epidemic in progress the basic reproduction number is typically estimated from 
incidence data and can vary significantly between different pandemic and seasonal outbreaks. 
Estimates for the 2009 HlNl pandemic are in the range 1.3 — 1.7 (Yang et al. 2009) and 
an average seasonal value is approximately 1.3 (Chowell et al. 2008). However, we have 
developed an explicit expression for the transmission rate per droplet (3d which allows a 
direct estimation of the reproduction number. Firstly, it is important to note that the 
infectious agent is not the droplet but the pathogens it carries. Therefore, the transmission 
rate per droplet (3d will depend the transmission rate per pathogen, (3d = (3pqdNp, where 
Qd is the probability of deposition in the human respiratory tract. The minimum infectious 
dose required to transmit the infection is thus implicitly incorporated into the transmission 
rate. The transmission rate per pathogen (3p is determined from the contact rate Cd of a 
susceptible with a droplet and the probability pd that such a contact will result in successful 
transmission (3p = CdPd- To derive the contact rate with a droplet it is assumed that each 
infected person is surrounded by a droplet cloud with volume Vd- It is further assumed that a 
susceptible individual comes in contact with a droplet through breathing during an encounter 
with this droplet cloud. If the average breathing rate is B and Tct is a characteristic time 
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parameter value 



c 


contact rate 


13 per day 


Pp 


pathogen concentration in the lung fluid 


3.71 X 10^ pathogens cm~^ 


B 


breathing rate 


24 m'^ per day 


Vcl 


personal-cloud volume of an infected person 


8 m^ 


Vd 


infection probability by an inhaled pathogen 


0.052 


Tct 


characteristic breathing (contact) time 


20 min 




pre-evaporation (spherical) droplet volume 


2.68 X 10-1° cm3 


qd 


inhaled droplet deposition probability 


0.88 



Table 1: Parameter values used to determine the transmission rate per droplet. 

of breathing during the encounter then the contact rate q can be expressed as q = Cy-^Td, 
where c is the average number of total contacts a susceptible individual has per unit time. 
The transmission rate per droplet is thus 

/3d = c—TctPdqdNp, 

and the number of pathogens per droplet can be determined by Np = VdPp, where Vd is the 
volume of the (spherical) pre-evaporative droplet and Pp is the pathogen concentration of 
the lung fluid. All the relevant parameter values for influenza are discussed by Stilianakis 
and Drossinos (2010), and the corresponding values for a 4/im droplet are summarized in 
Table [H The transmission rate per pathogen can now be calculated as 0.028/day and the 
transmission rate per droplet then evaluates to (3d = 2.45 x lO^^/day. Finally, the basic 
reproduction number can be estimated as Rq ~ 1.3, in agreement with seasonal estimates 
from the literature. Using the parameter values discussed above the dimensionless numbers 
are estimated as 

Ro = 1.3, p ^ 0.005, ^ ^ 0.0002, ^ 0.003. (20) 

The lifespan of droplets (and pathogens) is much less than that of an infected person 
such that p ^ 1. Also, for influenza, <^ ■?/' ^ 1 and the repopulation of the susceptible 
class following an outbreak is a slow process, primarily driven by loss of immunity and not 
natural births. The numerical solution approaches the endemic equilibrium with damped 
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Figure 5: Numerical solution of ( TT7|) -( IT^ . with initial conditions S{0) = 0.99, /(O) = 0.01 
and -D(O) = 10~^. Parameter values for influenza are Rq = 1.3, p ~ 0.005, (p ~ 0.0002 and 
ijj ^ 0.003. 

oscillations of large period ~ 10 years. Figure [S](a). The magnitude of the period is highly 
sensitive to the duration of immunity, with much shorter times observed between outbreaks 
if 1/cr is reduced. The small-time behavior of the model, where droplets rapidly achieve a 
balance with the infected population, is shown in Figure [5](b). After approximately 6 hours 
droplets and infected individuals achieve a balance and droplet dynamics thereafter follow 
the infected population. 



The model can also be used to describe the waterborne route of transmission (Codego, 
2001, Tien and Earn, 2010). An example of a pathogen that causes diarrhoeal disease 
in humans is the protozoa Cryptosporidium. Large outbreaks of Cryptosporidium due to 
contaminated drinking water are common (Eisenberg et al. 1998, Glaberman et al. 2002). 
The variable P{t) will represent the number of oocysts in a water source. This pathogen has 
been recorded surviving in an aqueous suspension for up to 12 months in cold temperatures 
a = 0.003/day, (Peeters et al. 1989). An infection typically persists for 6 — 9 days, fi = 
0.14/day, with infected individuals shedding up to 10^"^ oocysts per gram of faeces (Medema 
et al. 2009). With an average human faecal production rate of 106 g/day (Cummings et 
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al. 1992), we set k = 10 ''/day. The magnitude of the transmission rate by the waterborne 
route (in the absence of direct person-to-person transmission) has been estimated to be in 
the range [10~^^, 10~^] and we take /3 = 10~^°/day (Eisenberg et ah 1998). No immunity is 
acquired following infection and, for convenience, we assume that immunity is lost one day 
after recovery a = 1/ day. Using these parameter values we estimate dimensionless numbers 
as 

i?o = 2.38, p = 4:6.7, = 2.7xlO-^ ip = 7.1. 

The long lifespan of the Cryptosporidium oocysts yields p ^ 1. The solution of the model 
is plotted in Figure El An initial rapid period can be observed during which P is approxi- 
mately constant and S and / rapidly approach an intermediate state (5**,/*) = (0.97,0.02). 
The observed growth in the susceptible population in this region is an artifact of the 
choice of initial condition. At large times the solution converges to the endemic state 

Ee = (^e,/e,Pe) = (0.42,0.5,0.5). 

3 Seasonal variation in pathogen inactivation 

In this section we consider the dynamics of the model presented in Section [2] when seasonal 
forcing is applied to the pathogen inactivation rate. We make no assumption as to the pre- 
cise source of the forcing, merely that it manifests through increased pathogen survival in 
winter. We consider the specific case of the airborne transmission of respiratory infections 
by aerosol droplets, as described in Section 12.31 

3.1 The quasi-steady state 

To describe the seasonal cycle of respiratory infectious diseases we consider the effects of 
applying a seasonal forcing to the droplet model of Section 12.31 The forcing is implemented 



20 




Figure 6: Numerical solution of (l2])-(l4j), with initial condition S{0) = 0.9, /(O) = 0.1 and 
P(0) = 0.01. Parameter values for Cryptosporidium are Rq = 2.38, p ~ 46.7, ~ 2.7 x 10~^ 
and ~ 7.1. 



21 



via the droplet removal rate. The dimensional model is 



^^d(N-S) + a{N-S-I)-'^DS, 
— ^Kdl- a[t)D. 



We take 

a = 9 + + 7]i cos(cjt)}, 

where uj is the frequency of oscillation of the inactivation rate, 9 is the constant rate at which 
gravity and inhalation remove droplets. The parameter 770 is the pathogen inactivation rate 
in the absence of seasonal forcing and r]i represents the amplitude of the applied forcing 
term. For convenience, we denote = 9 + rjo. We nondimensionalise by scaling 

1 

S,I N, D -N, q; ~ ao, t 



ceo fJ> + a 

and we can write the dimensionless system as 



^ = ^(l-S) + ^{l-S-I)-Ilt,DS, (21) 

| = «,CS-/. (22) 

P§ = I-aD, (23) 

a = l + aiCos{nt), (24) 



where dimensionless numbers are defined as 



Ro = — 7 — —w, 0= — W = — -3, P= = 7fr> = > 



(25) 
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In the absence of seasonal forcing {ai = 0) the long-term behaviour is characterized by an 
endemic equilibrium state given by 

_ J_ , _ n _ (0 + V^)(i?o-l) 
'~Ro' ' ^ i?o(^ + l) ' 

in agreement with (jH]) where a = Oq. This equilibrium exists and is stable when Rq > 1 and 
one of either (p or ip are nonzero, enabling the repopulation of the susceptible compartment. 

The lifespan of airborne respiratory pathogens is relatively small, as gravity and pathogen 
inactivation invariably remove the pathogen-loaded droplets on a timescale significantly 
shorter than the infectious period, ^ T/, and we can reasonably assume p <^ 1. We 
adopt the quasi-steady state assumption and approximate the number of droplets by 

/ 

D = -. 

a 

In addition, unless immunity is permanently acquired following infection, the duration of 
the immune period is small compared with the average human lifespan and (p <^ tp. We thus 
allow repopulation of the susceptible class solely through the loss of immunity and neglect 
the natural birth process, i.e. 0^0. The model reduces to two equations 

§=Hl-S-I)-^IS, (26) 
at a 

'^ = -IS-I, (27) 

at a 

with a = 1 + cos(f2t). In contrast to models with sinusoidal forcing applied directly to the 
transmission rate /3, we find that the forcing appears in the denominator of the nonlinear 
term. Singularities are possible if a = 0, however, a small forcing amplitude ?7i <^ 1 always 
yields ai ^ 1 and division by zero is avoided. 
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3.2 Small-amplitude seasonal forcing 

Equation ( 1271) allows 5* to be written in terms of / as 

«(/ + /) 



S 



where the dot denotes the time derivative. This can be substituted into ( 126|) to yield a single 
equation for /, 



all + 



al — ai + aipl + RqP 



i + P[a- tpRo + tpRoI + RqI + aip] = 0. 



Assuming the seasonal forcing is small (ai <^ 1), a perturbation $,{t) to the endemic equi- 
librium Je such that 

I ^ Ie + ai^it) = Isit), 
satisfies the second order linear inhomogeneous equation 

^ + (V- + RoQ^ + HRo- l)e = h sin (nt) - ^ COS (nt)] . (28) 
dr^ at 

The characteristic equation of the corresponding homogeneous problem is 

X^ + itP + RoQX + HRo - 1) = 0, 
which has two negative roots Ai,2 < given by 



Al 



2 



The solution of the inhomogeneous equation is then 

^m=r -lAil* ,r -lA.lf I J WeRo -n'- sin (Qt) - leRom^ + 1) + COS {Qt) 

(29) 



24 



where Ci and C2 are arbitrary constants. At large times, ^ {t) approaches a periodic solution 
with period VL and amplitude 



A = L 



In addition, the solution for ^ has a phase $ which will satisfy 

Thus, for small-amplitude seasonal forcing, the infected population can be approximated by 

I^{t) = I^ + aiAsi\i{VLt + ^). (31) 

Linear theory indicates that the infected population oscillates about the endemic state Je 
with an amplitude aiA. The period of oscillation is identical to that of the applied seasonal 
forcing and independent of the duration of acquired immunity. However, the disease incidence 
and seasonally forced pathogen inactivation are out of phase with each other. This is readily 
seen by rewriting ( l3Ti) in the form 

J,(t) = Je + «iAcos|fi (^t + ^^^^|, (32) 

and there is a time shift of (2$ — 7r)/2fi between the two quantities, whereby /^(t) achieves 
its peak value before a does when 2$ — tt > 0. Furthermore, the maximum of Is{t) will not 
occur before the minimum of a{t) when 

2$ - vr TT , , 

Of course, a peak in a corresponds to the point of maximum pathogen inactivation and thus 
to a minimum in pathogen survival, and logically peak disease incidence should not coincide 
with minimal survival. Therefore, a more interesting and informative quantity is the delay 
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between minimum a (peak pathogen survival) and peak disease incidence. We denote this 
quantity as 5 and find 

where $ can be determined from ( l30i) . Thus, the extent of the observed delay between peak 
survival and peak incidence is determined from the specific disease parameters and the pe- 
riod of the seasonally forced inactivation rate. 



Another consideration in the relative positions of the maximum incidence and inactivation 
is the time of pathogen introduction, corresponding to the time the first infected individual 
is introduced into the population. For example, pathogen introduction at the point of peak 
survival could have radically different dynamics to an introduction at the time of peak 
inactivation. To address the consequences of the introduction time we consider a seasonal 
variation of the form 

a{t) = 1 + «! cos [^{t-t^)], (35) 

and then investigate the dependence of our results on to- Equation (135|) implies that the 
infected initial condition /(t = 0) may be introduced at different points during the seasonal 
variation of the pathogen inactivation rate. Namely, the infected person(s) at t = may be 
introduced in times that correspond to different pathogen inactivation rates. Given equation 
( l35ll . the calculations were repeated, and we found that the time shift t — to appears in 
the inhomogeneous term of the differential equation (!28|) and the infected population then 
satisfies 

J^(t) = + ttiAsin [n{t - to) + $] , 

which can be written as 

/,(t) = /e + «iAcos|fi(t-to + ^^^)}- (36) 

It can then be easily shown that the delay S does not depend on to, i.e. equation (I51|) 
remains valid, even for to 7^ 0. Therefore, we conclude that, while the solution of the 
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differential equation will invariably depend on the time of pathogen introduction, the delay 
between incidence I(t) and inactivation a(t) does not. 

3.3 The case of influenza 

As discussed in Section 12.31 in the case of influenza we find and it is reasonable to 

neglect the natural birth process. The values of the dimensionless parameters Rq, p and 
ip are given by equation (!20!) . The dimensionless frequency and amplitude are defined in 
( l25l) . Assuming a seasonal forcing with a period of one year yields a frequency of = 0.086. 
Droplet removal through inhalation is negligible and 6 can be approximated by the gravita- 
tional settling rate alone, for a 4/im droplet this yields 6 = 28.8/day (Robinson et al. 2012). 
The pathogen inactivation rate is assumed to be independent of droplet size and is taken to 
be riQ = 8.64/day (Hemmes et al. 1960). Finally, assuming a value of t]i = 0.1 for the forcing 
amplitude yields ai = 0.023. In reality, the value of rji depends on the seasonal driver (e.g. 
temperature, humidity) and its related seasonal cycle. This work merely assumes that rji 
is small and the chosen value determines the amplitude of the forcing term and thus the 
amplitude of the disease incidence. It is clear from ( 13T1) that the delay and period of the 
linearised solution are not affected by the choice of rji. 

The disease incidence Is(t) will reach its maximum before the maximal viral inactivation 
a{t) provided 2$ — tt > 0. From fl30|) . we calculate tan$ = 0.0167, the inverse of which 
yields $/ = 0.0167, where the subscript / indicates that this angle lies in the fundamental 
interval [— 7r/2, 7r/2]. However, the trigonometric coefficients in ( 129|) are both negative and $ 
must lie in the third quadrant, yielding $ = 0.0167 + vr. It follows that 2$ — vr = 3.175 > 0. 
Furthermore, the relative positions of the incidence and inactivation curves can be further 
restricted, by requiring the satisfaction of f l33|) . to ensure that the maximum of Is{t) will not 
occur before the minimum of a{t), and this condition reduces to 
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Figure 7: Influenza seasonality with initial conditions S{0) — 0.99, /(O) ~ 0.01 and -D(O) = 
10-3. Parameter values are Ro = 1.3, p ^ 0.005, (f) ^ 0.0002, ip ^ 0.003, Q = 0.086 and 
ai = 0.023, with to = 0. 
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Figure 8: Influenza seasonality with initial conditions S{0) = 0.99, /(O) = 0.01 and -D(O) = 
10-3. Parameter values are Ro = 1.3, p ^ 0.005, ^ 0.0002, tp ^ 0.003, Q = 0.086 and 
CKi = 0.023, with to = n/Q. 
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which is always satisfied by $/. 

A numerical solution of the dimensionless system ( l2T|) -f l2^ . showing the pathogen inacti- 
vation a{t), the disease incidence I{t) and the solution of the linearised system Is{t), is shown 
in Figure [71 The maximum disease incidence clearly occurs before maximum inactivation 
and after minimum inactivation. The time delay between peak virus survival (minimal a) 
and peak disease incidence is clearly visible and can be calculated from (IMj) as 5 ~ 18.07. 
Both curves, I(t) and a{t), are observed to oscillate with the same period, as predicted by 
the linear analysis. The hnearized quasi-steady solution Is(t), given by (|36|) . provides a very 
good approximation to the solution of the full system. Figure 7 displays the disease inci- 
dence for the case to = 0, such that the initial infected individual was introduced into the 
population at the point of maximum viral inactivation (minimum survival). To demonstrate 
the independence of the delay 6 on the time of viral introduction, the simulation was re-run 
with the infected individual introduced at the point of minimum inactivation (maximum 
survival), i.e. to = Tr/fi, Figure 8. The timing of peak disease incidence changes, however 
the predicted delay remains constant. 

A numerical solution is shown in terms of dimensional variables in Figure [9l where it is 
assumed that the peak in disease incidence occurs in January and, for illustration purposes, 
we have fixed the total human population at 6 million. The solution indicates that disease 
incidence peaks approximately three months (~ 90.33 days calculated from (IMI) ) after the 
influenza virus inactivation reaches its minimum, at which time virus survival is at its peak. 
The magnitude of this delay is dependent on the disease parameters. The delay 6 is plotted 
in dimensional form for varying Rq in Figure [TOT b). As Rq increases the delay decreases 
and the disease incidence peak occurs earlier after the peak in influenza virus survival. The 
amplitude of the infectious wave as a function of Rq is plotted in Figure [TOTa). Clearly, 
amplitude increases with increasing Rq, indicating the occurrence of more severe outbreaks 
for large _Ro values. 
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Figure 9: Numerical solution obtained with initial condition S'(O) = 0.99, /(O) = 0.01 and 
D{0) = 10-3. Parameter values Rq = 1.3, p ^ 0.005, ^ 0.0002, ^ ^ 0.003, fl = 0.086 and 
CKi = 0.023. All variables are dimensional. 






Figure 10: (a) Dimensional amplitude (= NaiA) as a function of Rq (b) Dimensional delay 
(= + d)) as a function of Rq. Parameter values = 6 x 10®, d = 4 x 10-^/day, 

/X = 0.2/day, 9 = 28.8/day rjo = 8.64/day, r/i = 0.1, oo = 
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4 Discussion 



In this work we have presented a model for the transmission of an infectious disease by 
means of a pathogenic state capable of surviving in an environmental reservoir outside of its 
host organism. For a reproduction number exceeding unity, i?o > 1, the pathogen dynamics 
are found to depend on the magnitude of a dimensionless number p, which represents the 
ratio between the pathogen lifespan and that of an infected individual. If the duration of 
the infectious disease state is significantly longer than the lifespan of the pathogen (e.g. for 
influenza), yielding p ^ 1, a rapid transient period is observed during which pathogen and 
infected numbers rapidly achieve a balance. Thereafter, the pathogen dynamics identically 
follow those of the infected population. At large times, the solution approaches a stable en- 
demic disease state, exhibiting damped oscillations as it does so. In contrast, if the pathogen 
is capable of surviving significantly longer than the duration of the infectious disease state 
(e.g. Cryptosporidium), yielding p ^ 1, the susceptible and infected populations first ap- 
proach an intermediate state during which the pathogen number remains approximately 
constant. However, at large times all variables converge to a stable endemic state. 

The effect of incorporating seasonal forcing through the pathogen inactivation rate was 
considered for the specific case of a respiratory infection mediated by airborne pathogen- 
loaded droplets. A quasi-steady state approximation allows the droplet dynamics to be 
neglected resulting in a classic SIRS model with the seasonal forcing appearing in the de- 
nominator of the transmission term. The assumption of small- amplitude seasonal forcing 
yields a linearized approximation to the infected population, given by f l32|) . It was found 
that the infected population oscillates about the endemic disease state with an identical 
period to that of the seasonally forced pathogen inactivation rate and the period is inde- 
pendent of the duration of acquired immunity. In addition, the peak in pathogen survival 
(corresponding to minimal pathogen inactivation) was found to precede peak disease inci- 
dence. For the particular case of influenza with a reproduction number of Rq = 1.3, this 
phase shift was calculated to be approximately three months, with shorter delays observed 
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for larger Rq values, Figure [TOl 

A final topic worth noting in relation to the seasonal model presented in Section [3] is 
the interpretation of the basic reproduction number. In the seasonal context, this no longer 
corresponds to the average number of secondary infections per primary infection since this 
value will invariably be dependent on time. Attempts have been made to define Rq in this 
context (Bacaer and Ait Dads 2012, Grassley and Fraser 2006) and the various attempts to 
derive a time-dependent expression have focused on the analysis of a time- dependent trans- 
mission term. However, in this work we attempt to tackle the issue of seasonality from a 
novel perspective by linking the seasonal signal to a palpable physical process. As such, the 
time-dependent Rq would be intrinsically linked to the pathogen inactivation rate which is a 
complex function of environmental factors, which in turn are highly variable with geographic 
location. To this end, an elaboration on the issue of time-dependence is beyond the scope of 
this work. 

The exact mechanism driving the seasonal occurrence of infectious diseases is uncertain. 
Several factors undoubtedly impact the dynamics, including human behavior and pathogen 
survival, and many of these factors are difficult to quantify. In particular, the burden placed 
on health systems by seasonal influenza, and the constant threat of deadly pandemics, makes 
understanding the transmission mechanisms of the virus an important issue. Experimental 
studies have identified the impact of environmental conditions on virus survival as a crucial 
factor. However, the primary sources of viral inactivation (temperature, humidity, solar 
radiation, etc.) have been difficult to identify, and indeed a complex combination of these 
sources working in unison is a likely culprit. More rigorous experimental work is needed to 
ascertain the primary source of the observed seasonality. Our model provides an insight into 
how annual variations in virus inactivation can influence disease incidence and the model 
is sufficiently general that it can be applied to various transmission pathways. A better 
understanding of pathogen removal processes could provide the opportunity to develop more 
complex models incorporating the effect and interaction of multiple environmental factors. 

32 



References 



[I] Anderson, R.M., May, R.M., 1981. The population dynamics of microparasites and their 
invertebrate hosts. Phil. Trans. Roy. Soc. London B. 291, 451-524. 

[2] Anderson, R.M., May, R.M., 1991. Infectious diseases in humans: dynamics and control. 
Oxford University Press. 

[3] Atkinson, M.P., Wein, L.M., 2008. Quantifying the routes of transmission for pandemic 
influenza. Bull. Math. Biol. 70, 820-867. 

[4] Bacaer, N., Ait Dads, E.H., 2012. On the biological interpretation of a definition for the 
parameter Rq in periodic population models. J. Math. Biol. 65, 601-621. 

[5] Bertuzzo, E., Casagrandi, R., Gatto, M., Rodriguez-Iturbe, I., Rinaldo, A., 2010. On 
spatially exphcit models of cholera epidemics, J. R. Soc. Interface. 7, 321333. 

[6] Burie, J.B., Calonnec, A., Ducrot, A., 2006. Singular perturbation analysis of travelling 
waves for a model in phytopathology. Math. Model. Nat. Phenom. 1, 49-62. 

[7] Cannell, J.J., Zasloff, M., Garland, C.F., Scragg, R., Giovannucci, E., 2008. On the 
epidemiology of influenza. Virology Journal. 5, 29. 

[8] Chowell, G., Miller, M.A., Viboud, C., 2008. Seasonal influenza in the United States, 
Prance and Austraha: transmission and prospects for control. Epidemiol. Infect. 136, 
852-864. 

[9] Codego, C.T., 2001. Endemic and epidemic dynamics of cholera: the role of the aquatic 
reservoir. BMC Infect. Dis. 1, 1. 

[10] Cummings, J.H., Bingham, S.A., Heaton, K.W., Eastwood, M.A., 1992. Fecal weight, 
colon cancer risk, and dietary intake of nonstarch polysaccharides (dietary flber). Gas- 
troenterology. 103, 1783-1789. 

[II] Dowell, S.F., 2001. Seasonal variation in host susceptibihty and cycles of certain infec- 
tious diseases. Emerg. Infect. Dis. 7, 369-374. 

33 



[12] Dushoff, J., Plotkin, J.B., Levin, S.A., Earn, D.J.D., 2004. Dynamical resonance can 
account for seasonality of influenza epidemics. PNAS. 101, 16915-16916. 

[13] Eisenberg, J.N., Seto, E.Y., Colford Jr, J.M., Olivieri, A., Spear, R.C., 1998. An analysis 
of the Milwaukee Cryptosporidiosis outbreak based on a dynamic model of the infection 
process. Epidemiology. 9, 255-263. 

[14] Fisman, D.N., 2007. Seasonality of infectious diseases. Annu. Rev. Public Health. 28, 
127-143. 

[15] Glaberman, S., Moore, J.E., Lowery, C.J., Chalmers, R.M., Sulaiman, I., Elwin, K., 
Rooney, P.J., Millar, B.C., Dooley, J.S., Lai, A.A., Xiao, L., 2002. Three drinking-water 
associated Cryptosporidiosis outbreaks. Northern Ireland. Emerg. Infect. Dis. 8, 631-633. 

[16] Crassly, N.C., Eraser, C, 2006. Seasonal infectious disease epidemiology. Proc. Roy. 
Soc. B. 273, 2541-2550. 

[17] Harper, C.J., 1961. Airborne micro-organisms: survival tests with four viruses. J. Hyg. 
(Lond). 59, 479-486. 

[18] Hemmes, J.H., Winkler, K.C., Kool, S.M., 1960. Virus survival as a seasonal factor in 
Influenza and Poliomyelitis. Nature. 188, 430-431. 

[19] Joh, R.I., Wang, H., Weiss, H., Weitz, J.S., 2009. Dynamics of indirectly transmitted 
infectious diseasees with immunological threshold. Bull. Math. Biol. 71, 845-862. 

[20] Keeling, M.J., Rohani, P., Crenfell, B.T., 2001. Seasonally forced disease dynamics 
explored as switching between attractors. Physica D. 148, 317-335. 

[21] Li, S., Eisenberg, J.N.S., Spicknall, I.H., Koopman, J.S., 2009. Dynamics and control of 
infections transmitted from person to person through the environment. Am. J. Epidemiol. 
170, 257-265. 

[22] Lofgren, E., Fefferman, N.H., Naumov, Y.N., Corski, J., Naumova, E.N., 2007. Influenza 
seasonality: underlying causes and modeling theories. J. Virol. 81, 5429-5436. 

34 



[23] Lowen, A.C., Mubareka, A., Steel, J., Palese, P., 2007. Influenza virus transmission is 
dependent on relative humidity and temperature. PLoS Pathogens. 3, 1470-1476. 

[24] Medema, G., Tennis, P., Blokker, M., Deere, D., Charles, P., Loret, J.F., 2009. Risk 
assessment of Cryptosporidium in drinking water. World Health Organization. 

[25] Miller, M.W., Thompson Hobbs, N., Tavener, S.J., 2006. Dynamics of prion disease 
transmission in mule deer. Ecological Applications. 16, 2208-2214. 

[26] Nelson, R.J., Demas, G.E., 1996. Seasonal changes in immune function. Quart. Rev. 
Biol. 71, 511-548. 

[27] Nicas, M., Nazaroff, W.W., Hubbard, A., 2005. Towards understanding the risk of 
secondary airborne infection: emission of respirable pathogens. J. Occup. Environ. Hyg. 
2, 143-154. 

[28] Nieuwhof, G.J., Conington, J., Bishop, S.C., 2009. A genetic epidemiological model to 
describe resistance to an endemic bacterial disease in livestock: application to footrot in 
sheep. Genetics Selection Evolution. 41, 19. 

[29] Peeters, J.E., Mazas, E.A., Masschelein, W.J., Villacorta Martiez de Maturana, I., 
Debackcr, E., 1989. Effect of disinfection of drinking water with ozone or chlorine dioxide 
on the survival of Cryptosporidium parvum oocysts. Appl. Environ. Microbiol. 55, 1519- 
1522. 

[30] Reluga, T., 2004. A two-phase epidemic driven by diffusion. J. Theor. Biol. 229, 249-261. 

[31] Robinson, M., Stilianakis, N.I., Drossinos, Y., 2012. Spatial dynamics of airborne infec- 
tious diseases. J. Theor. Biol. 297, 116-126. 

[32] Schulman, J.L., Kilbourne, E.D., 1963. Experimental transmission of Influenza virus 
infection in mice II: some factors affecting the incidence of transmitted infection. J. Exp. 
Med. 118, 267-275. 



35 



[33] Shaman, J., Kohn, M., 2009. Absolute humidity modulates influenza survival, trans- 
mission, and seasonality. PNAS. 106, 3243-3248. 

[34] Shaman, J., Pitzer, V.E., Viboud, C, Grenfell, B.T., Lipsitch, M., 2010. Absolute 
humidity and the seasonal onset of influenza in the continential United States. PLoS 
Biology. 8, el0000316. 

[35] Stilianakis, N.I., Drossinos, Y., 2010. Dynamics of infectious disease transmission by 
inhalable respiratory droplets. J. Roy. Soc. Interface. 50, 1355-1366. 

[36] Stohr, K., 2002. Influenza - WHO cares. Lancet Infect. Dis. 2, 517. 

[37] Stone, L., Olinky, R., Huppert, A., 2007. Seasonal dynamics of recurrent epidemics. 
Nature. 446, 533-536. 

[38] Tien, J.H., Earn, D.J.D., 2010. Multiple transmission pathways and disease dynamics 
in a waterborne pathogen model. Bull. Math. Biol. 72, 1506-1533. 

[39] Truscott, J., Fraser, C, Cauchemez, S., Meeyai, A., Kinsley, W., Donnelly, C.A., Ghani, 
A., Ferguson, N., 2012. Essential epidemiological mechanisms underpinning the transmis- 
sion dynamics of seasonal influenza. J. R. Soc. Interface. 67, 304-312. 

[40] Weber, T.P., Stilianakis, N.I., 2008. Inactivation of influenza A viruses in the environ- 
ment and modes of transmission: a critical review. J. Infect. 57, 361-373. 

[41] Yang, Y., Sugimoto, J.D., Halloran, M.E., Basta, N.E., Chao, D.L., Matrajt, L., Potter, 
G., Kenah, E., Longini Jr, l.M, 2009. The transmissibility and control of pandemic 
influenza A (HlNl) virus. Science. 326, 729-733. 



36 



